METHOD FOR IMAGE DESCRIPTION USING COLOR 



AND LOCAL SPATIAL INFORMATION 




BACKGROUND OF THE INVENTION 



5 



The present invention relates to a method for describing an image 



based on the color content of the image. 

Image description is a process for describing an image based upon 
the outcomes of the application of preselected measures to the image. 
Image description is useful in a number of applications such as digital 

10 image libraries where the descriptors are used as a basis for image 

indexing and retrieval. For image description to be practical and effective 
the outcome of the application of the measures to the image should be: 
(1) sufficient to distinguish between different images, (2) invariant to 
certain types of transformations of the image, (3) insensitive to noise, (4) 

15 easy to compute and (5) compact. Various methods of image description 
have been used and proposed with resulting image descriptors exhibiting 
these attributes to differing degrees. 

A paper by Swain et al. entitled COLOR INDEXING describes the 
use of color histograms to describe images. A color histogram of an 

20 image is obtained by calculating the frequency distribution of picture 
elements or pixels as a function of pixel color. Color histograms are 
invariant to translation or rotation of the image about the viewing axis. 
Color histograms can differ markedly for images with differing features. 
However, all spatial information about the features in the image is 

25 discarded in the creation of the color histogram. Therefore as long as two 
images have the same number of picture elements of each color it is not 
possible to distinguish between them using color histograms. This is true 
even if the two images contain features of completely different size or 
shape. For example, the total areas of the like colored (like hatched) 

30 geometric features of the two images of FIG. 1A and FIG. 1B are equal 
and require the same number of picture elements. The images cannot be 
distinguished on the basis of their color histograms even though the 
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features are clearly very different in size and number, and the images are 
easily distinguishable by the human eye. 

Several methods have been proposed to improve different aspects 
of the performance of color histograms. Strieker et al. in the paper entitled 
5 SIMILARITY OF COLOR IMAGES proposed the use of color moments. 
Color moments are statistical measures of the shape and position of the 
population distribution of pixel colors. In particular the color moments 
include a mean, a standard deviation and a skewness. Expressing the 
information contained in the color histogram in terms of a color moment 

10 results in a very compact image descriptor. Funt et al. in the paper 
entitled COLOR CONSTANT COLOR INDEXING proposed using the 
ratios of color triples [the red, the green and the blue pixels (RGB)] from 
neighboring regions of an image to reduce the effects of intensity 
variations. Rubner et al. in the paper entitled NAVIGATING THROUGH A 

1 5 SPACE OF COLOR IMAGES proposed the use of color signatures which 
is a plot of clusters of similar colors in an RGB color space. Using color 
signatures reduces the amount of data necessary to describe an image 
compared to that required for a color histogram. These methods improve 
some aspects of the performance of the image descriptors over the color 

20 histogram. However, like the color histogram, no spatial information is 
preserved. 

Several processes have been proposed which attempt to preserve 
some of the spatial information that is discarded in the construction of a 
color histogram. Pass et.al in the paper entitled HISTOGRAM 

25 REFINEMENT FOR CONTENT BASED IMAGE RETRIEVAL proposed 
refining the color histogram with color coherence vectors. In this process 
the coherence of the color of a picture element in relation to that of other 
picture elements in a contiguous region is determined. Even though the 
number of picture elements of each color is equal and, therefore, the color 

30 histograms are identical for two images, differences between features in 
the images will mean that the numbers of picture elements of each color 
which are color coherent will vary. Color coherence vectors do embed 
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some spatial information in the descriptors. Unfortunately, they require at 
least twice as much additional storage space as a traditional histogram. 

Rickman et al. in the paper entitled CONTENT-BASED IMAGE 
RETRIEVAL USING COLOUR TUPLE HISTOGRAMS proposed image 
5 description by construction of a histogram of the color hue at the vertices 
of randomly located triangular color tuples. Since the vertices of the 
triangular tuples are spaced apart, some spatial information is retained. 
Unfortunately, it is difficult to determine the dominant color of an image 
from the color tuple data. Further, the retained spatial information is 

10 difficult to interpret in a normal sense, therefore making it difficult to use 
the information for indexing an image database. 

"Color correlograms" were proposed for image description by 
Huang et al. in the paper entitled IMAGE INDEXING USING COLOR 
CORRELOGRAMS. A color correlogram quantifies the probability that a 

1 5 pixel of a particular color will lie at a specified radial distance from a pixel 
of a particular color in the image. The color correlogram provides a 
technique of measuring color coherence at different scales or distances 
from a point on the image. However, it is difficult to determine the 
dominant color of the image from a correlogram and it is difficult to 

20 interpret the correlogram in any usual human sense. 

Smith et al. in the paper entitled QUERYING BY COLOR REGIONS 
USING THE VISUALSEEK CONTENT-BASED VISUAL QUERY SYSTEM 
describe a method of image description using regions of color. Color data 
is transformed and the colors of the image are quantized and then filtered 

25 to emphasize prominent color regions. "Color set" values are extracted 
and a histogram is approximated by retaining those color set values above 
a threshold level. This method of image description 
requires image segmentation, a process that is difficult and 
computationally intensive. The region representation is rigid and variant to 

30 rotation or translation of images. 

"Blobworld" is a method of image representation proposed by 
Carson et al. in the paper entitled REGION-BASED IMAGE QUERYING. 
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In this method the image is segmented into a set of localized coherent 
regions of color and texture, known as "blobs." The "blobworld" 
representation of the image is the result of recording the location, size, 
and color of the segmented color blobs. This method provides 
5 considerable spatial information about the image, but the "blobworld" 
representation is rigid and variant to rotation or translation of images. 
Further, the image segmentation process is difficult and requires 
substantial computational resources. 



10 BRIEF DESCRIPTION OF THE DRAWINGS 

FIGS. 1A and 1B illustrate two images with features of different 

scale. 

FIG. 2 illustrates an image containing features of different colors or 
textures with delineated test areas. 
1 5 FIG. 3 illustrates the image of FIG. 2 with test areas of larger scale 

delineated on the image. 

FIG. 4 is an image for description with four square and four circular 
features. 

FIG. 5 is an image for description with a single square feature and 
20 a single circular feature where each feature has an area equal to the four 
features of the same geometric shape in FIG. 4. 

FIGS. 6A and 6B illustrate two similar images having features of the 
same size and shape but which have been translated and rotated. 

FIG. 7 is an exemplary illustration of the resulting image data for a 
25 first aspect of the present invention. 

FIG. 8 is an exemplary illustration of the resulting image data for a 
second aspect of the present invention. 

FIG. 9 is a graph of a nonbinary thresholding technique. 
FIG. 10 is an exemplary illustration of the resulting image data for a 
30 third aspect of the present invention. 

FIG. 11 is an exemplary color structure histogram. 
FIG. 12A illustrates an image with highly coherent color. 
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FIG. 12B illustrates an image with highly incoherent color. 
FIG. 13A illustrates an image with an 8x8 structuring element at 
single spacing. 

FIG. 13B illustrates an image with an 8x8 structuring element at 
5 double spacing. 

FIG. 14A illustrates a color space with quantization A. 

FIG. 14B illustrates a color structure histogram of FIG. 14A. 

FIG. 14C illustrates a color space with quantization B. 

FIG. 14D illustrates a color structure histogram of FIG. 14C. 
10 FIG. 14E illustrates a color space with quantization C. 

FIG. 14F illustrates a color structure histogram of FIG. 14E. 

FIG. 15 illustrates an image with two iso-color planes, P and Q. 

FIG. 16 illustrates an image with a single iso-color plane, PQ. 

FIG. 17 shows an exemplary data structure for colorQuant. 
15 FIG. 18 shows a HMMD color space. 

FIG. 19 shows an exemplary selection of available color spaces. 

FIG. 20 shows an exemplary order of the color spaces of FIG. 19. 

FIG. 21 illustrates one example of bin unification. 

FIG. 22 illustrates a technique for re-quantization and comparison. 
20 FIG. 23 illustrates linear pixel count versus code values. 

FIG. 24 illustrates non-linear pixel count versus code values. 

FIG. 25 illustrates one exemplary implementation of a color 
structure histogram descriptor extraction process. 

FIG. 26 illustrates one exemplary comparison for a query and a 
25 database descriptor. 

FIG. 27 illustrates an exemplary HMMD color space quantification. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
In existing systems of image description, the color or texture is 
30 quantified for a plurality of areas of predefined size and shape. The areas 
are preferably located on the image according to a predefined plan. The 
color or textural data for these areas of the image or statistical data related 
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thereto obtained are useful in describing the image and in distinguishing 
between images. The data obtained from each image may be referred to 
as an image descriptor. 

FIG. 2 illustrates the general application of image description using 
5 a generalized color histogram descriptor (characterization) based on an 
image having a triangular feature 2 of a first color and a circular feature 4 
of a second color. The color of the remainder of the image is a third 
background color. A number of square test areas 6 have been delineated 
on the image. The size and shape of the test areas may correspond to 

10 the size and shape of a predefined spatial structuring element 

encompassing a plurality of picture elements or pixels. While the spatial 
structural element defining the test areas illustrated in FIG. 2 is a square, 
there are no restrictions on the shape or size of the element. Regular 
shapes such as rectangles or circles may be more convenient in many 

1 5 applications than an amorphous shape or "blob." Also, the test area may 
be a scattered pattern of picture elements or pixels, akin to a shotgun 
blast. Likewise, the plan for locating the test areas on the image is not 
restricted to the rectilinear pattern illustrated in FIG. 2. 

A number of the test areas 6 of FIG. 2 lie wholly within the 

20 triangular feature 2. The color of the image in these test areas is the 

homogenous first color. Likewise, a number of test areas lie wholly within 
the circular feature 4 or the background. Over these test areas the image 
color is homogenous and can be quantified as either the second color or 
the background color, respectively. To varying degrees the remaining test 

25 windows overlap two or more regions of color. The colors in these areas 
are not homogeneous. 

Like the shape of the test areas and the plan for locating test areas, 
the size of the test area may be modified. Spatial information about the 
image is embedded in the data or image descriptor because the test areas 

30 have scale, that is, the areas encompass a plurality of picture elements. 
As can be seen by comparing FIGS. 2 and 3 changing the scale of the test 
area changes the number of test areas of each color. 
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Likewise if the sizes of the individual color regions of two images 
differ, the number of test areas of each color will likely vary. For example, 
the total areas of the four square 10 and circular 12 features of the image 
of FIG. 4 are equal to those of the square 20 and circular 22 features of 
5 the image of FIG. 5. As a result, the distribution of the population of 
picture elements as a function of color would be identical for the two 
images. However, as a result of the differences in sizes of the individual 
color regions of the images the number of test areas of each 
homogeneous color varies when the scale of the test area is held 

1 0 constant. In FIG. 5 there are more test areas that are the color of the 
circular feature than the test areas of FIG. 4 that lie wholly within the 
circular features. An image containing large uniform color regions or 
"blobs" will produce more test areas with the homogeneous color of those 
blobs than an image with smaller more scattered regions of color. 

15 While some test areas may lie completely within a region of 

homogeneous color, several of the test areas of FIG. 2 overlap two or 
more color regions. As a result the colors in these test areas are not 
homogeneous and must be quantified in some way to be useful in 
describing the image. For example, the mean values of the individual red, 

20 green, and blue (RGB) pixels, a transform of the RGB pixel values, or the 
mean color or the vector sum of the RGB intensity values might be used to 
describe the color of a test area of heterogeneous color. Since each test 
area having a heterogeneous color is likely to overlap two or more color 
regions to a degree differing from that of any other test area, there are 

25 likely to be as many mean colors or combinations of pixel intensities as 
there are test areas of heterogeneous color. Mapping the possible input 
values into a smaller number of quantized levels may be used to reduce 
the number of colors. For example, the RGB color data might be 
represented as the population of test areas in which percentage 

30 contributions of the red, green, and blue colors lie within certain ranges. 

As can be seen in FIGS. 2 and 3, only a small number of test areas 
may fall completely within the bounds of an image feature and, therefore, 
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have truly homogenous color. However, in several cases (see FIG. 2) a 
substantial part (less than all) of a test area is a particular color. The 
number of test areas included in the set of areas with homogeneous color 
can be increased by including in the application a test of homogeneity that 
5 would include in the data areas of "substantial" homogeneity. Likewise, 
accepting regions of images which are substantially homogeneous may be 
necessary for images which do not include many homogeneous color 
regions. 

For example, a test of homogeneity can be based on the standard 
1 0 deviation of colors of the picture elements in the test area. If <J k is the 
standard deviation of the pixel values in color channel k within a test area 
E then homogeneity can be defined by: 

H(e)=1-Ew k a k 

1 5 where w k is the weight coefficient for color channel k. 

An alternative homogeneity test function can be based on principle 
component analysis. A matrix A is defined as A=(Pjj) MxN where p tj is the /th 
color component of the rth pixel within a test area s. The singular values 
20 of A are determined by singular value decomposition. Letting p^, where 
k = 1 ,2, denote the singular values of A in descending order of 
magnitude, then homogeneity can be defined as: 

H(e)=1-gw k p k /p, 

where w k is the weight coefficient corresponding to singular 
value p k , /c>1. 

30 Data produced by the application of the image description can be 

incorporated into statistical representations which are familiar in the field. 
A "color blob" histogram can be constructed to present the frequency 
distribution of the population of test areas as a function of their color. For 
a given image I, a color blob histogram is the population distribution of all 

35 test areas of scale S, where S is the size of the test area in picture 
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elements. The color blob histogram is defined as an array h s that has an 
element h s c for each quantified color c belonging to the set C, that is 
C e C, and: 

5 h 8iC =|{ee| 8 |c( C )=c}|, 

where C is the set of all quantified colors and l s is the set of 
all color blobs of size s in the image I. 

10 The population distribution of test areas as a function of color can 

also be described by color blob moments which are the statistical 
moments of the color blob histogram. The color blob moments are 
extremely compact image descriptors. For a given image I, the first, 
second, and third statistical moments of the population distribution of the 

1 5 test areas of size S in each color channel k are: 

the mean (p) (first moment): 
p s k = 1 c k (£) 

the standard deviation (a) (second moment): 

the skew (A) (third moment): 

A s>k = 1 E (c k (e) - Ms , k ) 3 ) 1 ' 3 

30 

where: c k (e) is the kth color component of C(£). 

Referring to FIG. 7, the data resulting from a processed image may 
be represented as a set of quantized colors, |J 0 - |J 10 , together with an 
35 indication of the number of test areas of sizes S having a sufficiently 

homogeneous color matching one of the quantized colors. In other words, 
if |J 5 is red and six test areas of size X1 are sufficiently homogeneously 
red then |J 5 , the entry for |J 5 and S=X1 , would have a total of six. The 
result is a histogram where each of the entries totals the number of test 



20 



25 
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areas of size X1 having sufficiently homogeneous colors, as opposed to 
the summation of the colors of the individual pixels. The image may be 
processed with different test area sizes, S, to provide additional data. The 
resulting data from many images may be used for image comparison 
5 purposes. 

Image description using spatial test areas may result in invariance 
to rotation or translation of image features. In the two images in FIG. 6 
the square feature is translated vertically and horizontally while the 
triangular feature is rotated ninety degrees. The number of test areas 
1 0 having the homogeneous color of each of these features is unchanged. It 
can be shown for isotropic color areas that the color blob histograms and 
color blob moments are invariant to translation or rotation of image 
J~ features. 

fU The system may describe images on the basis of their texture or 

15 surface appearance. While color is a point property and can be described 
£ 4 by color histograms or other representations of the color properties of 

5 picture elements, texture is a local neighborhood property and texture 

fU descriptors describe the properties of an area surrounding a picture 

*= element. The texture of the individual test areas can be expressed in 

y 20 terms of mean texture descriptors, such as anisotropy, orientation, and 

contrast. The texture descriptors can be statistically described by a 
texture blob histogram. For an image I, a texture blob histogram for test 
areas containing S picture elements is the population distribution of test 
areas of size S, defined as an array h s that has an element h s t for each 
25 quantized texture model t contained in T and 
h s , t = |{£<= l s |t(e) = t}| 

where T is the set containing all quantized texture models. 
For a given image I, the texture blob moments for test areas of 
scale S are the first, second, and third statistical moments of the frequency 
30 distribution of the test areas of size S in each texture band k, that is: 
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the mean (fj) (first moment); 




the standard deviation (a) (second moment): 



5 




A sk = ( _1_ E (t k (8)- Msk ) 3 ) 

l"sl £6,S 



10 



where t k (e) is the kth component of t(e). 



The aforementioned technique counts the total number of test 
areas that are sufficiently homogeneous based upon the standard 
deviation of the color or texture. Unfortunately, selection of the threshold 
value for the standard deviation is difficult. If the threshold value is zero 

15 then no test area will likely be sufficiently homogeneous. Alternatively, if 
the threshold value is large then many of the test areas will likely be not 
very homogeneous, yet still be counted. FIG. 8 illustrates the percentage 
color distribution for the quantized colors for each test area. The resulting 
matrix has the number of occurrences of each quantized color as a 

20 function of color and color percentage. It is noted that the 100 percent 
column in FIG. 8 is the same as a single column of the aforementioned 
technique shown in FIG. 7. 

Referring again to FIGS. 2-5, the description of the technique is 
illustrated for matters of convenience as a set of test areas spaced apart 

25 from one another. To increase the invariance to translation and rotation 
the technique may involve locating the test area in an overlapping fashion 
at each pixel within the image. 

The size of the test area can have a profound effect on the number 
of sufficiently homogeneous test areas. Referring to FIGS. 4 and 5, if the 

30 test area used was selected to be larger than the square and circular 

features 10 and 12 (FIG. 4) but less than the square and circular features 
20 and 22 (FIG. 5), then processing FIG. 4 may result in no sufficiently 
homogeneous regions. However, processing FIG. 5 would result in 
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several sufficiently homogeneous regions. In this manner the differences 
in the number of sufficiently homogenous test regions would be increased 
which would allow for easier differentiation between images using such 
measures. 

5 The technique described herein is applicable to any suitable color 

space, such as Y/Cb/Cr. The pattern and size of the test areas on the 
images may be changed or be random, if desired. 

The aforementioned homogeneity test provides a result that is 
either sufficiently homogenous (yes or "1") or not sufficiently homogenous 

1 0 (no or "0"), in a manner similar to a step function. Such a homogeneity 
test is sensitive to noise because slight variations in the standard 
deviation, which is a calculated quantity, may change the result of the 
homogeneity test if the standard deviation is close to the threshold. 
Accordingly, the aforementioned homogeneity test is sensitive to noise 

15 and doesn't take into account finer gradations in the amount of 
homogeneity. Referring to FIG. 9, the homogeneity threshold 
determination may include a "soft" thresholding mechanism. The 
thresholding mechanism provides a floating point measure (e.g., not a 
binary yes/no determination) of the homogeneity in reference to some 

20 measure of the homogeneity, such as the standard deviation. The 
thresholding mechanism may provide a gradual increase in the 
homogeneity as the standard deviation decreases. In this manner small 
changes in the standard deviation, in a region proximate the threshold, will 
not result in significant changes in the measure of the homogeneity. In 

25 addition, the particular selection of the threshold value is less critical to 
achieving accurate results. Other non-binary functional definitions of the 
homogeneity as a function of some measuring criteria may likewise be 
used, if desired. 

Referring again to FIG. 8, the percentage color distribution for the 
30 quantized colors for each test area is illustrated based on an equal 

percentage distribution for each column. However most images contain a 
large variety of color content in most regions of the image. Accordingly, 
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the color distributions for most images tend to be distributed toward the 
smaller percentages. In other words, typical images contain relatively few 
large regions of substantially pure homogenous color. With relatively few 
significant regions of homogenous color, the portions of the matrix with 
5 larger percentage values tend to be primarily zero which wastes space 
and does not provide an effective technique of discriminating between real 
world images that contain smaller differences. Referring to FIG. 10, to 
overcome these limitations and to maintain a relatively compact matrix, the 
matrix may include smaller percentage ranges at the smaller percentages, 
10 with increasing percentage ranges toward the larger percentages. This 
maintains a small matrix, which is suitable for embedded systems, while 
providing more accurate discrimination between images with similar color 
content. 

It is to be understood that the aforementioned description regarding 
1 5 a "soft" thresholding technique and modified matrix is likewise applicable 
for texture. 

The present inventors considered the aforementioned techniques 

and realized that the selection of the percentages, such as shown on 

FIGS. 8 and 10, is at least partially arbitrary. In addition to being arbitrary, 
♦ 

20 if the selection of quantized colors (p x ) are finely spaced, generally 
resulting in a large number of available quantized colors, then minor 
changes in the colors of the image as a result of noise will significantly 
change the overall result. Further, in addition to the arbitrary percentages 
and the effect of finely quantized colors (p x ), if the percentages are finely 

25 spaced then slight differences in the amounts of colors will result in 
substantial differences in the resulting image descriptors. As it may be 
observed, it becomes increasingly more difficult to accurately compare 
different image descriptors derived from different but visually similar 
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images because of the susceptibility to variances in the image descriptors 
themselves. 

In contrast to attempting to further refine the percentages and 
available quantized colors, the present inventors postulated that if the 
5 percentage boundaries, as shown in FIGS. 8 and 10, are eliminated and 
the homogeneity test is simplified so that it merely determines if any of the 
quantized colors exist within the test areas, then a robust and readily 
usable image descriptor is achievable. Referring to FIG. 1 1 , the indicies 
(e.g., 0-255) along the axis represent a quantized color in the chosen color 

10 space, thereby forming a color structure histogram. To create the color 
structure histogram each different color contained in the image within each 
test area (or a selected set thereof) is identified. Then each different 
identified color is quantized according to the quantized color regions. The 
duplicate quantized colors are discarded for each test area. In other 

15 words, each quantized color in the test region of the image for each test 
area is counted merely once. The resulting color structure histogram is a 
one-dimensional histogram with the data contained therein representing 
more than merely the total color distribution of the image. The additional 
information contained in the resulting color structure histogram includes, 

20 for example, the frequency of the colors and the color coherence of each 
quantized color (spatial information). In essence, the system de- 
emphasizes the effect of spatial regions of coherent color in the histogram 
and emphasizes the effect of spatially incoherent regions of colors in the 
histogram. Referring to FIGS. 12A and 12B, the color histogram can 
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distinguish between two images in which a given color is present in 
identical amounts but where the structure of the groups of pixels having 
the color is different in the two images. For example, FIG. 12A would 
record a value of 90 (9x10) for the color in the histogram. In contrast FIG. 
5 12B would record a value of 459 (nine for each interior color (9x45), three 
for each edge color (3x4), and six for each color one away from the edge 
(7x6)) for the color in the histogram. A comparison between a traditional 
color histogram and a particular implementation of the color structure 
histogram, illustrating the benefits, is described in ISO/IEC JTC 1/SC 
10 29/WG 11/M5572, Maui Hawaii, December 1999, incorporated by 
reference herein. 

A DDL representation syntax for the color structure may be defined 
as follows: 

<complexType name= M ColorStructureType n > 
1 5 <complexContent> 

<extension base=" VisualDType"> 
<sequence minOccurs— ' 1" maxOccurs— 'l f, > 
<element name— 'Values" minOccurs=" 1 " maxOccurs-' 1 "> 
<simpleType> 

20 <list itemType="unsigned8 M > 

<minLength value="3/32'7> 
<maxLength value="256'7> 
</list> 
</simpleType> 
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</element> 
</sequence> 

<attribute name="colorQuant" type= n mpeg7:unsigned3" 

use"required f 7> 
5 </extension> 

</complexContent> 
</complexType>. 

The retrieval effectiveness of the color structure histogram is 
1 0 significantly better than that of the traditional histogram, for descriptors of 
the same number of "bins" (i.e., number of quantized colors). The color 
structure histogram is particularly effective in comparison to the traditional 
histogram when descriptors with a small number of bins are compared, 
i.e., the case of coarse color quantization. 
1 5 The extraction complexity of the color structure histogram is as 

follows. If K is the number of quantized colors in the histogram, and S is 
the number of pixels in the structuring element, then the order of 
complexity is 0(S+K) per pixel, where 0() generally refers to the order of 
computational complexity operator, well known in the art as so-called big 
20 "O" or "Landau" notation. The complexity of computing the histogram over 
the entire image is 0((S+K)n), where n is the number of pixels in the 
image. Assuming color quantization is performed prior to histogram 
extraction, only integer summations, multiplications, comparisons, and 
memory read/writes are needed to compute the color structure histogram. 
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If the number of bins in the histogram is n, then the order of 
complexity of histogram matching is O(n), in cases when an ^ distance is 
used as a similarity measure, where ^ refers to an {<, norm (sum of the 
absolute differences). If the ^ distance is used, only integer summation, 
5 comparison operations, and memory read/writes are needed to match two 
color structure histograms. 

After further consideration of the test areas an attempt was made to 
determine the optimal size of a test region. It is to be understood that the 
optimal test size determination may likewise be used for other types of 

10 histograms that incorporate spatial information. It is problematic to 
determine an optimal test size with respect to retrieval accuracy for the 
structuring element. One of the difficulties is that a fixed size structuring 
element is not optimal for all images. After processing two different 
images representing the same scene at different scales using the same 

1 5 sized test area the present inventors were surprised to observe that the 
resulting color structure histograms, normalized to take account of the 
differing image sizes, were very different. This would not be the case with 
the traditional histogram. After observing this unexpected result, the 
present inventors then postulated that the primary source of the difference 

20 were the different scales of the two images. Based upon these 

postulations and observations, the present inventors then determined that 
the size of the test area (or equivalently the structuring element) should be 
modified in accordance with the size of the image being processed. 
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Accordingly, a relatively larger image should use a relatively larger test 
area, whereas a smaller image should use a relatively smaller test area. 

An analysis of a database of images with approximately the same 
size (e.g., 320x240 and 352x288) using structuring elements (test areas) 
5 of different sizes, different pixel densities, and different layout patterns of 
positions within the image was performed. The structuring elements used 
were 1x1, 2x2, 4x4, 8x8, and 16x16. The 1x1 structuring element is a 
special case which is equivalent to extracting a traditional color histogram. 
The test results suggest that the retrieval performance generally improves 

10 with increasing structuring element size (having a given pixel density and 
given layout pattern). Significant performance improvements may be 
observed when increasing the structuring element size from 1x1 (regular 
histogram) to 2x2, and to 4x4, and to 8x8. In many cases, the 
performance improvement becomes small when increasing the structuring 

1 5 element further. The sensitivity of the performance to the size of the 
structuring element is relatively low (i.e., there is no clear performance 
"peak" for a particular structuring element size). The exact structuring 
element size (within a few pixels) does not appear to be critical, with an 
8x8 structuring element appearing to be preferable. Improvement was 

20 observed when the structuring element was increased by factors of two. 
After consideration of the retrieval accuracy data resulting from the 
database analysis, the present inventors determined that it is not 
necessary to precisely relate the structuring element size to the image, but 
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rather it is sufficient to use factors of two which allows a straightforward 
logarithmic-exponential relationship and limits computational complexity. 

While any technique may be used to modify the relative size of the 
structuring element, the preferred technique is described below. Referring 
5 to FIGS. 13A and 13B, the spatial extent of the structuring element should 
depend on the image size; however, the number of samples in the 
structuring element may be maintained constant, by sub-sampling the 
image and structuring element at the same time. The number of samples 
in the structuring element is preferably maintained at 64, layed out in an 

10 8x8 pattern, and the distance between two samples in this pattern 

increases with increasing image sizes. This technique is equivalent to 
sub-sampling the image by a power of two and then using a structuring 
element of 8x8 pixels. That is, the technique may be interpreted, in one 
embodiment, as resizing the image to a fixed base size and always using 

15 the same densely packet 8x8 structuring element. The technique may be 
performed "in place" in software, that is, the sub-sampling of the image 
may be done implicitly by simply skipping samples during processing, 
while computing the color structure histogram. The sub-sampling factor 
and the spatial extent of the structuring element width and height can be 

20 computed at the same time as follows. Let E be the spatial extent of the 
structuring element size, i.e., the spatial extent is preferably ExE. Let K be 
the sub-sampling factor to be applied, where K=1 implies no sub- 
sampling, K=2 implies sub-sampling by 2 horizontally and vertically, etc. K 
and E are preferably computed as follows: 
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p=max{0,round(0.5*log 2 (width*height)-8)} where K=2 P and 

E=8*K 

For example, an image of size 320x240 using the formula above 
will yield K=1 and E=8, in which case the structuring element is simply 8x8 
5 with no sub-sampling performed as shown in FIG. 13A. An image with 
size 640x480 using the formula above will yield K=2 and E=16, in which 
case the spatial extent of the structuring element is 16x16 and sub- 
sampling is 2x2 as shown in FIG. 13B, which results in a structuring 
element with spatial extent of 8x8 on the sub-sampled image. Note, that 
10 images smaller than 256x256 are a special case in the sense that K=1 
and E=8 in all cases. This avoids up-sampling smaller images to a bigger 
size and at the same time performs sufficiently well. 

An implementation of the variable sized test area, illustrating the 
benefits, is described in ISO/I EC JTC 1/SC 29/WG 1 1/M5785, 
15 Noordwijkerhout, the Netherlands, March 2000, incorporated by reference 
herein. 

It is desirable to have available descriptors of different length, i.e., 
different numbers of "bins". As previously described, this corresponds to 
descriptor extraction in a color space that has been more coarsely or finely 
20 quantized. In general, a small descriptor corresponds to a more coarsely 
quantized color space. However, the color space may be quantized in any 
non-uniform manner, if desired. The different sized descriptors permits 
the particular system to select, at least in part, the storage requirements 
necessary for storing the color structure histograms. In addition, the 
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selection of the size of the descriptors of the color structure histogram 
permits the system, at least in part, to determine the system's complexity 
and computational requirements. For example, with a limited number of 
images and nearly unlimited available storage, then a descriptor with a 
5 relatively large number of bins may be desirable. Where there is an 

unusually large number of images with limited additional available storage 
and limited computational resources, then a descriptor with a relatively 
limited number of bins may be desirable. For embedded systems where 
storage space is severely limited, a descriptor with a severely limited 

10 number of bins may be desirable. The available descriptors may be 

selected as desired, such as for example, 256, 200, 175, 130, 96, 75, 32, 
and 12. It is to be understood that multiple descriptor sizes may be used 
with any image descriptor system, including but not limited to color 
structure histograms. 

15 FIGS. 14A-F, describe the relationship between the quantized color 

space and the associated bin-layout along the independent axis of the 
color (or color structure) histogram descriptor. It also describes the 
relationship between two histograms derived from two different color 
space quantizations. A two dimensional color space divided into a small 

20 number of disjoint subsets, each encompassing a contiguous region of 
space, is shown in FIG. 14A for illustrative purposes only. In practice the 
dimensionality of the color space may be higher, typically being three and 
its shape may be arbitrary. Also in practice the number of subsets may be 
larger or smaller, their shape may be arbitrary, and the portions of space 



-21 - 




they contain may be highly disconnected, even consisting of one or more 
disconnected (discrete) points. To facilitate the discussion, these disjoint 
color space subsets shall be called "cells" although, as just mentioned, 
their shape and form may be arbitrary. FIG. 14A shows a particular 
5 quantization of the displayed color space which shall be denoted as "A" 
type quantization. By numbering the cells from 0 to N-1 where N (here 
N=16) is the total number of cells, and then numbering with the same 
numerals the bins of an N bin histogram, shown in FIG. 14B, a bijective 
relationship is established between the histogram bins and the color space 

10 cells. That is, each bin corresponds to one and only one cell and, 
conversely, each cell corresponds to one and only one bin. The 
assignment of the N numbers to both the color space cells and the 
histogram bins is arbitrary but in practice an orderly scheme such as that 
shown in FIGS. 14A-F is used. The value in a particular bin, say the kth 

1 5 bin, of the color structure histogram is determined, as discussed earlier, by 
the number of positions of the structuring element within the image that 
contain a color which is located within the kth color space cell. For a 
traditional histogram, the value in the kth bin is the number of times a pixel 
having a color in the kth cell occurs within the image. 

20 FIG. 14C illustrates a re-quantization of the color space, which shall 

be denoted "B" type color space quantization. By re-quantization, it is 
meant that the color space is partitioned into a different set of cells, 
possibly but not necessarily a different number of cells. The independent 
axis of the histogram associated to FIG. 14C is shown in FIG. 14D. FIG. 
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14A and FIG. 14C illustrate the case where there is little relationship 
between the quantization cells of a space and the cells of its re- 
quantization. If one is given the histogram of FIG. 14B and wishes to 
convert it into the histogram of FIG. 14D, for reasons of interoperability 
5 and without reference to the associated image itself, then the following 
difficulty arises. How, precisely, can the values in the histogram bins of 
FIG. 14B be combined to obtain bin values for FIG. 14D? Because of the 
bijective relationship between bins and color space cells, this is equivalent 
to asking how to re-apportion the number of pixels whose colors lie in cells 
1 0 of the "A" quantization to the number of pixels that lie in the cells of the "B" 

fit 

!| quantization. The difficulty is illustrated by considering the cell of the "B" 

%j quantization whose index is 3. This cell contains portions of cells 4, 5, 7, 

s and 8 from the "A" quantization shown in FIG. 14C by the dashed 

boundaries. Thus some portion of the number of pixels having a color 
□ 15 lying in each of these "A M quantization cells should contribute to the value 

•pas:. 

in bin 3 of the histogram of FIG. 14D corresponding to the "B" 
quantization. But without reference to the original image pixels this 
apportionment is difficult to determine. 

The inventors conducted experiments to test various possible 
20 schemes by which to do this apportionment rationally. One idea was to 
apportion pixels having color in a given cell of "A" type quantization to a 
given cell of "B" type quantization in proportion to the area of the given "A" 
cell which overlaps the given "B" cell. Retrieval results from using this 
method to re-quantize descriptors were poor because the method does 
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not (and cannot) take into account where in the given "A" quantization cell 
the pixel colors were originally located. The inventors realized that only in 
the case where cells from the "A" quantization lie completely inside or 
completely outside a cell of the re-quantized space could such an 
5 apportionment be made. For in that case, all or none, respectively, of the 
pixels in the given "A" quantization cell would, ipso facto, lie in the given 
cell of the re-quantized space. 

FIG. 14E shows a color space re-quantization of the "A" 
quantization which has this property. This is denoted the "C" quantization 

10 of the color space. Observe that every "A" cell lies completely inside or 
outside of some "C" cell. Equivalently, every "C" cell boundary is an "A" 
cell boundary. With such a re-quantization of the color space the 
derivation of the "C" quantization histogram values from the "A" 
quantization histogram values may proceed. A preferred technique of 

15 derivation is to combine by addition the values of those "A" histogram bins 
which correspond to "A" cells that have been combined into "C" cells by 
re-quantization. FIG. 14F illustrates this for two "C" histogram bins, bin 0 
and bin 3. Bin 0 of FIG. 14F corresponds to cell index 0 in FIG. 14E. This 
cell is the (trivial) combination of just one "A" quantization cell from FIG. 

20 14A, namely the cell with index 4. Hence the value placed in bin 0 of FIG. 
14F is derived solely from the value found in bin 4 of the "A" histogram of 
FIG. 14B. 

As a non-trivial example, consider bin 3 of the "C" histogram of FIG. 
14F. This corresponds to "C" cell index 3 of FIG. 14E which encompasses 
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precisely cells 1,2,3, and 1 1 from the "A" color space quantization. 
Hence the values from the "A" histogram found in bins 1,2,3, and 1 1 are 
combined, preferably by addition, to obtain the derived value for bin 3 of 
the "C" histogram in FIG. 14F. 
5 An exemplary example of how this re-quantization may be 

accomplished is described below for purposes of illustration. Let A be the 
color space quantization of a histogram and B be the target re- 
quantization. Let l A be a given color bin index in the A histogram. In HSV 
(hue-saturation-value) color space, for example, re-quantization may be 

*~ 10 performed by first mapping l A to Hq A , Sq A , and, Vq A , the quantization 

indices of the three HSV color components for the A type quantization. 

% 4 The ^mapping is defined by inverting the map that takes individual 

quantized color indices and delivers a histogram bin index. Next, the 

f5 three color indices are de-quantized according to: H=(Hq A +0.5)/nHq A , 

P 15 where nHq A is the number of levels to which H was originally quantized in 

the A type and where H is a floating-point quantity. The same formula, 
with suitable changes, applies to S and V. Then l B is computed by re- 
quantizing H, S, and V, according to the quantization levels of the B type 
quantization and re-computing the histogram bin index, l B , from Hq B , Sq B , 
20 and Vq Q . This defines a map from l A to l B . The histogram amplitude index 
in l A is simply added to l B . It can be shown that this is equivalent to adding 
the histogram amplitudes at l A and l B . 

While re-quantization may be applied to color histograms and color 
structure histograms, the present inventors came to the startling 
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realization that this is not an optimal operation to perform when using the 
color structure histogram descriptors for image retrieval, as described 
below. In particular, this is not an optimal operation when the color 
structure histograms are extracted at different quantization levels and then 
5 subsequently re-quantized. The principal reason for this behavior is in the 
nature of the color structure histogram and is closely related to the 
reasons why color structure histograms normally out-perform the 
traditional histogram. Referring again to FIGS. 12A and 12B, they 
illustrate qualitatively, the behavior of the color structure histogram in the 

1 0 presence of two pathological but instructive types of color structures within 
an iso-color plane, the plane of pixels all having the same color. In FIG. 
12A pixels of the same color, call it color P, are clumped together in a 
rectangular "blob". For the sake of description this dumpiness may be 
referred to as coherence. The more coherent an iso-color plane is, the 

1 5 more likely it is that groups of pixels within the iso-plane will be found 
close together. Conversely, the more incoherent the iso-color plane, the 
more its pixels will tend to be far apart, where "far apart" is with respect to 
the dimensions of the structuring element used for the color structure 
histogram. 

20 The coherence of FIG. 12A, neglecting edge effects, contributes 

(8+2)x(7+2)=90 counts to the (un-normalized) color structure histogram 
bin, the P-bin, that corresponds to color P. This is because a pixel of color 
P will be found within the structuring element at 90 different positions of 
the structuring element. On the other hand, the count for FIG. 12B will be, 
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neglecting edge effects, (8x7)x(3x)=504, because each pixel now 
contributes 9 total counts to the color structure histogram bin. 

The corresponding traditional histogram will have 56 (un- 
normalized) counts in either case. Accordingly, the traditional histogram is 
5 blind to the incoherence of the color structure whereas the color structure 
histogram, in addition to measuring the amount of each color, is also 
sensitive to the incoherence within the iso-color plane. This additional 
information is the principal reason why the color structure histogram out- 
performs the traditional histogram. Likewise, the present inventors 

10 realized this is also principally why the color structure histogram can not 
be expected to perform well under re-quantization, as explained below. 

Referring to FIG. 15, let A again denote the initial color space 
quantization and B a coarser scalable re-quantization. A second color, Q, 
is introduced which has the following three properties: (i) its structure is 

1 5 also incoherent; (ii) its pixels are spatially near the pixels of color P; and 
(iii) its position in color space is near enough to color P that it will lie in the 
same quantization bin as P, cell PQ, when re-quantized. Color Q also 
contributes 504 counts to its respective improved histogram bin, the Q-bin. 
The corresponding traditional histogram again gets (8x7)=56 counts in its 

20 Q-bin. 

Presume, for purposes of illustration, that the color structure 
histogram and the traditional histogram are re-quantized. The P-bin and 
Q-bin become the new PQ-bin. For the traditional histogram the count in 
PQ-bin is 112, the sum of counts in the P-bin and Q-bin, because that is 
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how one does scalable re-quantization: a bin in the B quantization gets the 
contents of those bins in the A quantization that it contains. Notice that 
this is the same value that would be in the traditional histogram PQ-bin if 
the image had started out with B quantization. This is because a pixel in 
5 the B space has color PQ if and only if it had color P or color Q in the A 
quantized space. In other words, re-quantization for the traditional 
histogram is additive, (or, more properly, homomorphic) in the sense that 
combining two colors into one and then counting it is the same as 
individually counting the two colors and then adding the results. 

1 0 The behavior is quite different for the color structure histogram. 

When the color structure histogram is re-quantized, one adds the counts 
in all the bins that map to a given re-quantized bin just as with the 
traditional histogram. This is the best that one can do in the absence of 
knowledge of the structure of the associated iso-color plane. The result is 

15 1008 counts. However, if the image starts out in the B quantized color 
space a very different result occurs. This can be observed in FIG. 16, 
where different color pixels have now become the same color. It may be 
observed that the incoherence of the iso-color plane is reduced in relation 
to FIG. 12B. Therefore, one can expect to get a lower count in the PQ-bin 

20 of the color structure histogram than resulted when re-quantizing the color 
structure histogram itself because re-quantizing can not take into account 
the color structure. In fact, the count would be 736 for FIG. 16, were the 
descriptor extracted from the image quantized in the B type color 
quantized space, given a 3x3 structuring element. 

-28- 




As a result, re-quantized color structure histograms are not 
homomorphic. A color structure histogram extracted from a B quantized 
image is significantly different with respect to the norm, from one that is 
re-quantized from A to B. Testing of the re-quantization of the color 
5 structure and traditional histograms is described in ISO/I EC JTC 1/SC 
29/WG 11/M6018, Geneva, May 2000, incorporated by reference herein. 

One of the attribute names within the MPEG-7 DDL definition of the 
descriptor presented earlier is colorQuant which specifies the color space, 
the color quantization operating point, and determines the number of 

10 ColorStructure values used in the DDL representation syntax. Its 
semantics may be specified as illustrated in FIG. 17. The variable, 
colorQuant, may take on suitable values, for example, 001 , 010, 01 1 , and 
100. The values field contains the ColorStructure descriptor data which is 
organized in an M element array of 8-bit integer values, h(m) for me{0, 1 , 

15 M-1}. The number, M, of bins may be chosen from the set {256, 128, 
64, 32} of allowable operating points. The bins of an M-bin descriptor are 
associated bijectively to the M quantized colors, c 0 , c^ c 2 , c^, of the M- 
cell color space, which is defined later. The value of h(m) represents, in a 
non-linear manner to be described, the number of structuring elements in 

20 the image that contain one or more pixels with color c m . 

It is to be understood that any color space may be used, as 
desired. However, for purposes of completeness the preferred color 
space is referred to as "HMMD M . The HMMD color space is defined by a 
non-linear, reversible transformation from the RGB color space. There are 
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five distinct attributes (components) in the HMMD color space. The 
semantics of the five attributes are defined as follows: 
Hue: Hue; 
Max: max(R, G, B); 

indicates how much black color the image has, giving 
the flavor of shade or blackness; 
Min: min(R, G, B); 

indicates how much white color the image has, giving 
the flavor of tint or whiteness; 
Diff: Max-Min; 

indicates how much gray the image contains and how 
close to the pure coir, giving the flavor of tone or 
colorfulness; 
Sum: (Max+Min)/2; and 

simulates the brightness of the color. 
Referring to FIG. 18, the HMMD color space has a double cone 
appearance consisting of blackness, whiteness, colorfulness, and hue. A 
selection of available color spaces may be ordered in any desired 
sequence, such as the sequence shown in FIG. 19. The available color 
spaces may be further represented as a binary value, if desired, such as 
the binary representation shown in FIG. 20. 

Normally the image descriptors are extracted and compared in a 
common color space. It is considerably more difficult to compare image 
descriptors that are derived from different color spaces. 

-30- 




In light of the realization that it is not optimal to re-quantize color 
structure descriptors for comparison with one another, the present 
inventors determined that the color structure histogram should always be 
initially extracted from the image at the finest quantization granularity, 
5 such as 256 levels. Referring to FIG. 21 , after extraction at the finest 
quantization the descriptor may be re-quantized by combining appropriate 
bins, such as by simple summation. In this manner the other levels, such 
as 128, 64, and 32 may be determined in a consistent and uniform 
manner which is independent of the color coherence of the image. 
ft 10 Referring to FIG. 22, the database of color structure histograms are 

!J created by initially quantizing each image at the highest quantization level, 

*-4 such as 256 bins, at block 250. The quantized images as a result of block 

~ 250 are then re-quantized to the desired number of bins, such as 128, 64, 

r ~ 

rll 

U and 32, at block 252. The search query is initially quantized at the highest 

p 15 quantization level, such as 256 bins, at block 260. The quantized image 

as a result of block 260 is then quantized to the desired number of bins, , 
such as 128, 64, and 32, at block 262. The quantized images as a result 
of blocks 250 and 260 need not be re-quantized, if desired. Block 270 
determines if the query descriptor is quantized at a different level than the 
20 particular histogram descriptor. If the two descriptors have the same 

number of bins then the descriptors are compared, at block 272. If the two 
descriptors have a different number of bins then the descriptor is re- 
quantized to match the quantization of the other descriptor, at block 274 
prior to comparison. The descriptors may both be re-quantized to the 
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same number of bins, if desired. The re-quantized descriptors, now being 
at the same size, are compared at block 272. With each color structure 
histogram being quantized to the same size, in the same manner, the 
color structure histograms will be consistent with one another and 
5 accordingly the spatial information contained therein will be uniformly 
treated. 

After further consideration of a histogram including spatial 
information, especially when each quantized color is merely counted once 
for each test area, a significant number of the bins contain relatively small 

10 numbers. To further reduce the storage requirements for the histogram, 
the bin amplitudes are quantized into a selected set of code values. For a 
color structure histogram the maximum value that any particular bin 
amplitude may obtain is a predefined number, namely, (N-S x +1)x(M-S y +1), 
where N is the horizontal width of the structuring element in pixels, M is 

15 the vertical height of the structuring element in pixels, S x is the horizontal 
width of the structuring element in pixels, and S y is the vertical height of 
the structuring element in pixels. It is noted that this maximum value is the 
same as the traditional color histogram, where S x = S y =1 . With the 
maximum potential value being known, the resulting histogram may be 

20 normalized in a well defined manner. Referring to FIG. 23, an example of 
an inter-relationship between the normalized total pixel count and the 
resulting code values is shown. Traditionally, the pixel count is uniformly 
quantized which includes a linear relationship between code values and 
quantized amplitudes, as shown by the diagonal dotted line in FIG. 24. 
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Referring to FIG. 24, another example of an inter-relationship between the 
normalized pixel count having a non-linear relationship to code values. 
This is an example of non-uniform quantization. 

Most of the data within typical color structure histograms are small 
5 numbers plus a few large numbers, such as illustrated by FIG. 1 1 . When 
comparing two histograms comprised mostly of small numbers, typically 
by the absolute difference of one histogram from another, the result will 
primarily be smaller numbers. This decrease in the apparent differences 
between the small numbers is further decreased by subsequent amplitude 
10 quantization, if performed. Accordingly, the remaining few large numbers 
will tend to dominate the comparison between two color structure 
histograms. To compensate for the tendency of large code values (i.e., 
large numbers) to dominate while small code values (i.e., small numbers) 
become nearly irrelevant, the present inventors determined that the 
15 amplitudes should be non-uniformly quantized which induces a non-linear 
relationship between amplitudes and code values. An exemplary 
distribution of the different code values may divide the bin amplitude range 
into six regions, and subsequently allocate a different number of 
quantization levels uniformly within each region. The thresholds to divide 
20 the bin amplitude range (between 0.0 and 1 .0) into 6 regions are (or 
approximately): 

ThO 0.000000001 ; (or a number significantly smaller than 

0.037 or substantially equal to zero) 
Th1 0.037; 
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Th2 0.080; 
Th3 0.195; and 
Th4 0.320. 



The number of quantization levels (or code values) allocated to 
5 each region are (or approximately): 



N0^ 


=1 


between 0.0 and ThO; 


Nl 


=25 


between ThO and Th1; 


N2: 


=20 


between Th1 and Th2; 


N3; 


=35 


between Th2 and Th3; 


N4: 


=35 


between Th3 and Th4; and 


H5 


=140 


between Th4 and 1 .0. 



The threshold values may be modified, as desired. 
In contrast to the traditional wisdom of uniformly quantizing the bin 
amplitudes, the improved technique uses a non-uniform amplitude 
15 quantization technique. An implementation of the non-uniform 

quantization of amplitudes is described in ISO/IEC JTC 1/SC 29/WG 
1 1/M5218, Beijing, July 2000, incorporated by reference herein. 

Referring to FIG. 25, one exemplary implementation of a color 
structure histogram descriptor extraction process is shown. A "raw" 256- 
20 bin histogram is accumulated (e.g., compiled) directly from the image, at 
block 300. At this point, bin amplitudes are un-quantized and reside in the 
"linear" domain, i.e., linearly related to the number of structuring elements 
that contain the color associated with the bin. If 256 bins are desired then 
block 302 branches to block 304 which non-uniformly quantizes the 
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amplitude of the bins, such as using the techniques previously described. 
If less than 256 bins are desired then block 302 branches to block 306 
which re-quantizes the color space by bin unification, such as using the 
techniques previously described. The result of bin unification at block 306 
5 is still in the "linear" domain. The results of block 306 are clipped at block 
308 to a maximum amplitude, which avoids integer "rollover" if a limited 
number of bits are used. The result of clipping by block 308 is provided to 
block 304 which non-uniformly quantizes the amplitude of the bins. The 
result of block 304 non-uniform amplitude quantization, which provides 

1 0 code values which are non-linearly related to the number of structuring 
elements that contain the color associated with the bin. After considerable 
analysis, the present inventors determined that the re-quantization via bin 
unification in the "linear" domain provides increased retrieval performance 
over bin unification in the "non-linear" domain using code values. This 

15 increased performance, the present inventors determined, is primarily the 
result of decreased clipping. 

Referring to FIG. 26, when a query and a database descriptor are 
presented for comparison to a similarity measure their sizes must agree. 
Given a database descriptor of size M 320 and a query descriptor of size 

20 N 322, the larger of the two descriptors is reduced in size to match the 
smaller of the two. The code values of the descriptor to be reduced are 
first converted to (quantized) linear amplitudes at block 326. The 
conversion of code values to linear amplitudes normally have the following 
properties: (i) there is a linear relationship between the resultant 
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amplitudes and the mid-interval values of the non-uniform quantization 
intervals within [0,1] defined previously, and (ii) these linear amplitude mid- 
interval values are represented by B bits, where B is preferably 20. The 
bin unification is performed at block 328. In particular, if it is assumed that 
5 M>N, then the mapping of the bins in the M-bin descriptor to the bins in 
the N-bin descriptor is defined by re-quantizing the color represented by 
each bin of the M-bin descriptor into the N-cell color space, and then 
computing the bin index that represents each re-quantized color. The 
result of block 326 is a descriptor with non-uniform amplitude quantization. 
10 During bin unification the sum of two bins are preferably clipped at block 
330 to the maximum possible linear amplitude, 2 B -1 . Then, the linear 
amplitudes of the reduced descriptor are converted back to non-linear 
code values. 

FIG. 27 shows a slice of the HMMD space in the diffsum plane for 
15 zero hue angle and depicts the quantization cells for the 128-cell operating 
point. Cut-points defining the subspaces are indicated in the figure by 
vertical lines in the color plane. The diff-ax\s values that determine the cut- 
points are shown in black at the top of the dashed cut-point markers along 
the upper edge of the plane. Horizontal lines within each subspace depict 
20 the quantization along the sum-axis. The quantization of hue angle is 
indicated by the gray rotation arrows around each cut-point marker. The 
gray number to the right of a rotation angle corresponds to the number of 
levels to which hue has been quantized in the subspace to the right of the 
cut-point. For example, Figure 14 states that the hue values associated 
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with the subspace between cliff = 60 and diff- 110 (i.e. subspace 3) are 
quantized to 8 levels. This agrees with the entry in Table 15. 

The bijective mapping between color-space cells and descriptor bin 
indices is given explicitly by the numbers within the cells. The ordering of 
5 these numbers is first from bottom to top (parallel to the sum-axis), then 
from cliff-sum plane to diff-sum plane (around the ftue-axis) staying within 
a subspace, and finally from subspace to subspace. For example, the 
cells of Figure 14 closest to the bottom edge in subspaces 2 and 3 are 
numbered 32 and 64. The jump is due to the fact that there are four sum 

10 levels and 8 hue levels for this subspace. The numbers within the 
subspace, therefore, increase from 32 to 32 + 4*8 - 1 = 63. 

The terms and expressions that have been employed in the 
foregoing specification are used as terms of description and not of 
limitation, and there is no intention, in the use of such terms and 

1 5 expressions, of excluding equivalents of the features shown and described 
or portions thereof, it being recognized that the scope of the invention is 
defined and limited only by the claims that follow. 
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